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Abstract. We present a comparative study of X-ray emission from a sample of active M dwarfs with spectral 
types M3.5 - M4.5 using XMM-Newton observations of two single stars, AD Leonis and EV Lacertae, and two 
unresolved binary systems, AT Microscopii and EQ Pegasi. The light curves reveal frequent flaring during all 
four observations. We perform a uniform spectral analysis and determine plasma temperatures, abundances and 
emission measures in different states of activity. Applying multi-temperature models with variable abundances 
separately to data obtained with the EPIC and RGS detectors we are able to investigate the consistency of the 
results obtained by the different instruments onboard XMM-Newton. We find that the X-ray properties of the 
sample M dwarfs are very similar, with the coronal abundances of all sample stars following a trend of increasing 
abundance with increasing first ionization potential, the inverse FIP effect. The overall metallicities are below 
solar photospheric ones but appear consistent with the measured photospheric abundances of M dwarfs like these. 
A significant increase in the prominence of the hotter plasma components is observed during flares while the cool 
plasma component is only marginally affected by flaring, pointing to different coronal structures. AT Mic, probably 
a young pre-main-sequence system, has the highest X-ray luminosity and exhibits also the hottest corona. While 
results of EQ Peg and EV Lac are presented here for the first time, AT Mic and AD Leo have been investigated 
before with different analysis approaches, allowing a comparison of the results. 
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1. Introduction 

Observations with ROSAT have shown the forma- 
tion of X-ray emitting coronae around late-type cool 
dwarf stars with outer convection zones to be universal 
llSchmitt et all 1 19951 ISchmitt & LiefkeL l2004h . Coronal 
structures on the Sun are dominated by their magnetic 
properties with closed structures, the coronal loops or 
loop arcades, containing predominantly the X-ray emit- 
ting, confined hot high-density plasma and being the lo- 
cation of flare events. In the Sun the interaction between 
the radiative and the outer convective zones powers a 
dynamo, leading to coronal activity. The ubiquitous oc- 
currence of X-ray emission among cool stars and their 
dependence on rotation strongly suggests also a mag- 
netic character of their activity. Observed activity lev- 
els of cool stars, as measured by the ratio of bolomet- 
ric and X-ray luminosities, span a range of values from 
~ 10~ 8 to ~ 10~ 3 that are correlated with dynamo ef- 
ficiency. The efficiency of a solar-like magnetic dynamo 
is characterised by the inverse of the so-called 'Rossby- 
Number', defined as the ratio between rotational period 
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and convective turnover time. This activity-rotation re- 
lation saturates at £x/£boi~ 10~ 3 , a level > 1000 times 
higher than that of the Sun. The similarity of stellar and 
solar X-ray emission and flares suggests similar basic phys- 
ical mechanisms, with more active stars showing more fre- 
quent and larger flares. Flaring has been proposed as a 
possible ma jor coronal he ating mechanism for the Sun via 
nanoflares jParkerl Il988|) as well as for other cool stars 
including M dwarfs. For recent reviews on this topic, see 
iKashvap et all l|2002|) or lGudel et ail (|2003|) . 

Active M dwarfs, belonging to the so-called flare stars, 
turned out to be particularly strong coronal X-ray sources. 
Frequent flaring on these objects was first discovered in the 
optical regime, later in EUV and X-rays. These late-type 
dwarf stars are generally more active than solar-like stars 
since they have larger spin-down times. They also typically 
show higher coronal temperatures. M dwarfs of spectral 
types later than approximately M3 - M5 are expected to 
be fully convective and therefore lack an interface between 
radiative and convective zones. The observed high levels 
of activity then require the onset of alternative dynamo 
mechanisms different to that of the Sun. 

The stars chosen here for analysis are known high X- 
ray emitters with spectral types in the range M3.5- M4.5 
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and are located in the immediate solar vicinity (distance 
5-10 pc). EQ Peg (Gl 896A, Gl 896B) is a visual bi- 
nary consisting of two M dwarfs of spectral type M3.5 
and M4.5, both of which are flare stars, as well as AT Mic 
(Gl 799A, Gl 799B), an active binary system consisting 
of a M4.5 and a M4 star. Both binaries are separated by 
roughly 30 AU, making interaction between the compo- 
nents quite unlikely. On the other hand, AD Leo (Gl 388) 
is a single M3.5 star similar to EV Lac (Gl 873). Our sam- 
ple M dwarfs are also all fast rotators ( y sin i > 5 km/s ) 
l|Delfosse et"allll998tlFuhrmeister et~al"l l2004|) in the sat- 
urated regime, i.e Lx/^boi ~ 10~ 3 . Stellar parameters im- 
portant for this work are summarized in Tabled Note that 
the spectral classification in the literature is not unique for 
all M dwarfs in our sample and may differ by one spec- 
tral subclass; the values listed here are taken from the 
SIMBAD database. 



bv lBesselaar et all l|2003|) and again an inverse FIP effect 
was found. An analysis of spatially resolv ed X-ray emis- 
sion o f the EQ Peg system was presented bv lRobrade et alJ 
( 2004). A comparative study on coronal density diagnos- 
tics with He-like tripl ets that includes w ork on our sample 
stars can be found in lNess et alJ ll2002[) and with speci al 
focus on the X-ray sizes of coronae in llSfess et alJ l)2004|h 
The above listed M dwarfs were observed with XMM- 
Newton and the goal of this paper is to present a compara- 
tive analysis of these X-ray observations and to determine 
those properties that are typical for active M dwarfs as a 
class and those attributed to individual objects. In Sect.[5] 
we describe the observations and the methods used for 
data analysis. In Sect.|3|we present the results subdivided 
into different physical topics which is followed by a sum- 
mary and our conclusions in Sect. 01 



Table 1. Stellar parame ters; vsini fro m Dclfo sse et alJ 
l|l998|) ex cept AT Mic llFuhrmeister et all l2004|) and 
EQ Peg llStauffer fc HartmannL 119861) . J from 2MASS 
poin t source catalog, Lb ni calcu lated from J using values 
of iKenvon fc Hartmannl l|l995l) . Lboi/Lx with Lx (this 
work, PN). 



Star 




d(pc) 


vsini (km/s) 


J (mag) 


logL bo! 


logL(, o! / Lx 


AT Mic 


A+B 


10.2 


11.7(A) 


5.81 


32.46 


-3.00 


EQ Peg 


A+B 


6.2 


12/18(A/B) 


5.79 


32.04 


-3.22 


AD Leo 




4.7 


6.2 


5.45 


31.94 


-3.21 


EV Lac 




5.0 


6.9 


6.11 


31.73 


-2.99 



While X-ray emission from AD Leo and AT Mic 
was already detected by HEAO 1 l|Kahn et all Il979h 
and EQ Peg i s conta ined in the Einstein Stellar Survey 
( Vaianaet all Il98llh EV Lac w as first detected by 
EXOSAT l|Schmitt & RossoLll988|) . All stars in our sam- 
ple were also detec ted in the ROSAT all-sky survey 
( Hii nsch et all Il999(l and were observed by various other 
X-ray missions. A comparative analysis of X-ray spectra 
from AD Leo and EV Lac obtain ed with BeppoSAX an d 
ROSAT PSPC was performed bv lSciortino et alJ |I1999Il 
iFavata et ail l)2000a|) presented a study of a n extreme X- 
ray flar e on EV Lac observed with ASCA. IFavata et alJ 
" 2000bh studied AD Leo data taken by Einstein, ROSAT 



and ASCA and found the coronal structures to be compact 
and the quiescent coronal luminosity remarkably constant 
over decades. 

Parts of the XMM-Newton data of AT Mic were previ- 
ously analysed with special focus on elemental abundances 
and emission measure s during flaring and quiescent state 
( Raassen et all l2003fl . While the coronae of low activ- 
ity stars like the Sun IjLaming et a,llll995h show an en- 
hancement of low-FIP elem ents, the so-cal l ed FI P (First 
Ionization Potential) effect, iRaassen et alJ ll2003h discov- 
ered in the corona of the highly active star AT Mic an 
inverse FIP effect, which appears to flatten during the 
flare state. A comparative study of flares on AD Leo ob- 
served with XMM-Newton and Chandra was performed 



2. Observation and data analysis 

The four targets, AD Leo, AT Mic, EV Lac and EQ Peg, 
were observed with XMM-Newton using varying detector 
setups and exposure times in the range of 15-35ksec. No 
attempts were made to separate the binaries in this work, 
therefore data from AT Mic and EQ Peg is a superpo- 
sition of the respective system components. A detailed 
description of the observations is provided in Table |2j 
Useful data were collected in all X-ray detectors on- 
board XMM-Newton, respectively the EPIC (European 
Photon Imaging Camera) and RGS (Reflection Grating 
Spectrometer) detectors which are operated simultane- 
ously. The EPIC instrument consists of three CCD cam- 
eras with two different CCD designs, respective two MOS 
and one PN, providing imaging and spectroscopy in the 
energy range from 0.15 to 15keV with good angular and 
moderate spectral resolution. In this work EPIC data is 
used only when the detectors observed in the imaging 
mode and therefore data from both MOS detectors are 
presented only for EQ Peg . The RGS is a grating spec- 
trometer with high spectral resolution consisting of two 
identical spectrometers. The RGS operates in the energy 
range from 0.35 to 2.5 keV (5 - 35 A) with a spectral resolu- 
tion of ~4.0eV (0.06 A, FWHM) at 1.0 keV in first order. 
A detailed description of th e XMM-Newton instruments 
is given bv lEhle et all |2f 



Table 2. Observation log 



Target 



Instrument (Mode) 



Du 



Obs. Time 



EQ Peg MOS (FF, thick F.) 14600 2000-07-09Tll:39-15:43 

PN (FF, thick F.) 12410 2000-07-09T12:20-15:47 

RGS (HER/SES) 15610 2000-07-09Tll:31-15:51 



AT Mic MOS1 (SW.med. F.) 25400 2000-10-16TOO:28-07:32 
PN (SW,med. F.) 25100 2000-10-16T00:42-07:40 
RGS (HER/SES) 28310 2000-10-16T00:20-08:ll 



AD Leo MOS2 (LW, thick F.) 35800 2001-05-14T20:55-06:52 
PN (SW,med. F.) 35000 2001-05-14T21: 11-06:54 

RGS (HER/SES) 36400 2001-05-14T20:48-06:55 



EV Lac MOS2 (SW.med. F.) 32260 2001-06-03T07:55-16:53 
PN (SW,med. F.) ' 31470 2001-06-03T08:ll-16:55 
RGS (HER/SES) 32860 2001-06-03T07:49-16:56 
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The data were reduced with the standard XMM- 
Newton Science Analysis System (SAS) software, version 
5.4.1 with calibration files updated in January 2004. Light 
curves and spectra were produced with standard SAS tools 
and standard selection criteria were applied for filtering 
Ehle et all 120041). Spe ctral analysis was done 
with XSPEC VI 1.3 llArnaudl ll99r3) . 

Data from the detectors of different type were fitted 
separately to allow comparison of the results and to check 
for possible cross-calibration problems. Spectral analysis 
of EPIC data is performed in the energy band between 
0.2- 12.0 keV. RGS first order spectra in the full energy 
range, i.e. 0.35- 2.5 keV (5 -35 A) are used. Data of the 
same detector type, e.g. RGS1 and RGS2, and in the 
EQ Peg case MOS1 and MOS2, were analyzed simulta- 
neously but not co-added. For spectral analysis the data 
was cleaned of proton flares by removing the affected time 
periods, and in the EQ Peg case the data was addition- 
ally cleaned for pile-up by excluding the inner part of the 
source emission. The applied selection criteria were ad- 
justed according to observation mode. The background 
was taken from source-free regions on the detectors. 

We apply the same type of spectral model to all data 
to perform a self-consistent analysis. Therefore we use the 
global fitting approach, where the whole spectrum is fitted 
instead of a line-base d method which wou ld be applica- 
ble only to RGS data. lAudard et"aTI l|2004j) applied both 
methods to data of YY Mensae taken by Chandra and 
XMM-Newton and found consistent results. These authors 
fitted data from RGS and MOS in selected energy ranges 
simultaneously, similar to our joint fits, which are briefly 
discussed in Sect 13.2. 21 Larger deviations were found when 
using data of different detectors, a problem also aff ecting 
our analysis. However, ISanz-Forcada et al. | ll2003h anal- 
ysed Chandra and XMM-Newton data of AB Dor, a highly 
active near-main-sequence K star, and found the line- 
based method to be better suited to fit their high reso- 
lution spectra, leading to slightly higher abundances and 
lower emission measures to be present in the derived mod- 
els com pared to the global fitting approach. iGiidel et all 
l|200lh analysed EPIC and RGS spectra from different 
observations of AB Dor using the global approach and 
found different absolute abundance values but a pattern 
(increasing depletion of abundance of low FIP elements up 
to iron followed by the inverse pattern for elements with 
higher FIP) overall consis tent with the results derived by 
ISanz-Forcada et all l(2003l ). demonstrating the robustness 
of abundance ratios in general. 

For the analysis of the X-ray spectra we specifically 
use multi-temperature models with variable but tied abun- 
dances, i.e. the same abundance pattern in all tempera- 
ture components. Such models assume the emission spec- 
trum of a collisionally-ionized diffuse optically-thin gas a s 
calculated with the APEC code fsee lSmith et all l200l]) . 
Abundances are calculated relative to solar photospheric 
values as given bv lAnders fc Grevessd ill 9891). For iron and 
oxyge n we use the updated values of iGrevesse fc Sauvall 
( 199§§. Comparing stellar coronal with solar photospheric 



abundances is of course crude, but precise photospheric 
abundances are unknown for the observed targets; a de- 
tailed discussion of the abundances is given in Sect . 13. 2. 
Due to the proximity of the stars, absorption in the in- 
terstellar medium is negligible at the wavelengths of in- 
terest and was not applied in our modelling. Our fit pro- 
cedure is based on x 2 minimization, therefore spectra are 
rebinned to satisfy the statistical demand of a minimum 
value of counts per spectral bin; 25 for data with an av- 
erage count rate above 5ctss _1 and 15 for data with an 
average count rate below 5ctss _1 and for all data used 
in time resolved analysis. Bad channels are always ex- 
cluded from the fit; these include channels at high en- 
ergy in the EPIC spectra that contain insufficient source 
counts (see Fig.^J). The first model uses three tempera- 
ture components (3 -T model) with free temperatures and 
emission measures. Models with additional temperature 
components were checked, but did not improve the fit 
results significantly. In order to approximate physically 
more realistic continuous temperature distributions of the 
plasma, we use a 6-T model on a logarithmic, almost 
equidistant grid (0.2,0.3,0.6, 1.2, 2.4, 4.8 keV), which cov- 
ers roughly the temperature range from 2 to 60 MK and 
samples those spectral regions where the XMM-Newton 
detectors are most sensitive. 

Consider next the different spectral resolution, AE/E, 
and energy range of the EPIC and RGS instruments. 
Especially at low energies, AE/E is quite large for the 
EPIC detectors, while in the RGS most of the strongest 
emission lines are separated in this energy range, thus per- 
mitting a more accurate determination of plasma param- 
eters. Some elements, in particular C and N, show promi- 
nent features only in the RGS. Therefore, to fit the EPIC 
data the abundances of these elements were set to the RGS 
value in order to avoid unphysical solutions. For elements 
such as Al, Ca and Ni, whose abundances in the solar 
photosphere are an order of magnitude lower than that of 
Fe and which show only very weak features in the X-ray 
spectrum, we adopted a coronal abundance of 0.4 times 
the solar photospheric abundance (a typical value for Fe 
in this sample). The EPIC detectors are able to measure 
higher energy X-rays, important for a reliable determina- 
tion of the properties of the hotter plasma components, 
with the MOS detectors providing a slightly better spec- 
tral resolution and the PN detector providing greater sen- 
sitivity. 

An important feature of the considered spectral mod- 
els is the interdependency between the emission measure 
(EM= J n e nijdV) of temperature components on the one 
hand and elemental abundances with temperature sensi- 
tive emission lines in the considered temperature range 
on the other hand. This is especially true for EPIC data 
with their moderate spectral resolution. Another general 
source of uncertainties is the incomplete knowledge of 
atomic data used to calculate the spectral models, see e.g. 
IPorauet et all l|200lh and references therein. 
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3. Results 

3.1. Investigation of light curves 

The PN light curves of the four observations are shown in 
Fig-ffl these light curves were extracted from a circle with 
50 " radius around each source and the temporal binning 
is 100 s for all sources. Light curves were cleaned (AD Leo, 
AT Mic, EV Lac) for obvious data dropouts and high pro- 
ton flare background was subtracted from the light curve 
(EV Lac). As the PN exposure of AD Leo began during 
a flare, we included scaled data from the MOS2 exposure, 
which started ~ 1 ksec earlier, to show the peak of the 
flare. 

The light curves of the individual targets show obvi- 
ous variability in all sources. The largest flares that oc- 
curred in these observations involve X-ray flux increases 
of factors 2-3, which is considered moderate in view of 
flux increases of factors 10-300 that have been observed 



flares on EQ Peg l|Katsova et all 12002ft and EV Lac 
l|Favata et all l2000aft . For an analysis of a larger flare on 
the M dwar f Prox ima Centauri with XMM-Newton see 
iGiidel et all (12004ft . While our observations indicate that 
flaring with flux increases of factors 2-3 appears regularly 
on timescales of a few hours, these larger events are much 
rarer. 

Weaker flaring (flux changes of around ten percent or 
lower) on timescales of a few minutes is visible throughout 
all four observations. In this phase the stars are in a state 
of low but permanent activity. This continual weak flaring 
makes it, in addition to the larger flares, difficult to deter- 
mine a minimum or baseline flux level. This flux or activ- 
ity level is denoted here as the quasi-quiescent level. It has 
been suggested that this emission can be interpreted as a 
superposition of a large nu mber of even smaller flares with 
a power law distr i bution ( Hudson! Il99lt iKashvao et all 
l2002Hr T iidel et a,lll2003ft . 

The flare light curves exhibit a variety of shapes and 
often show substructure. The observed rise times are only 
somewhat shorter than the decay times for the larger im- 
pulsive flares during the AT Mic and the EV Lac ob- 
servation, pointing to relatively compact flaring struc- 
tures. The different types of flares a s obser ved on the Sun 
are described in e.g. iFeldman et al 1 111995ft and references 
therein. The rise times from quasi-quiescent to flare max- 
imum are typically in the range of 1 ksec for the largest 
flares with the rise steepening towards maximum. The bi- 
naries are not resolved in these light curves, but all com- 
ponents are expected to be flar e stars. For EQ Peg it was 
shown bv lRobrade et al.l ^2004) that both stars contribute 
to the flaring, with EQ Peg A being brighter than B by 
a factor of 3-4. The separation of AT Mic is too small to 
resolve the components. 

Flaring involves significant amounts of energy release 
and matter transfer into the corona and may therefore in- 
volve changes in the physical properties of the coronal 
plasma, e.g. temperatures, emission measures, densities 
and abundances, a well known phenomenon from previous 



analyses of stellar flares. Therefore we introduce 'flaring' 
and 'quasi-quiescent' phases for each star; an interval is 
selected as flaring if the count rate of the time bins exceed 
the median count rate and the peak count rate within the 
interval exceeds a 'minimum flare value', chosen to sepa- 
rate clear flare events and flickering. The selection of two 
phases of activity is of course somewhat arbitrary, and 
observations with stronger flares may be separated into 
three or more phases of in- or decreasing activity. With 
our data this was only done for the AT Mic observation, 
which is separated in three states of activity, the quasi- 
quiescent phase (QQ), the medium activity phase and the 
outstanding large flare. The count rates, selection criteria 
and selected flaring periods are summarized in Table OH 



Table 3. Selection of flaring data (PN). 





EQ Peg 


AT Mic 


AD Leo 


EV Lac 


Median rate 


6.42 


9.40 


8.60 


7.58 


Min. flare rate 


7.0 


12.0 (high) 


9.5 


8.5 






10.0 (med) 






Flare intervals 


5.5-11.6 


14.7-19.1 (h) 


0.0-7.7 


0.5-2.3 


(ksec after 




9.7-12.1 (m) 


12.2-14.9 


8.4-12.2 


start of 




19.2-25.0 (m) 


24.1-25.7 


19.4-24.0 


observation) 






27.4-28.4 










32.0-33.8 





3.2. Spectral analysis 

3.2.1. RGS and EPIC spectra 

The high resolution RGS spectra of our sources in the 
wavelength range 5-35A are shown in Fig.|2j In this en- 
ergy range the X-ray spectra of our M dwarfs are remark- 
ably similar, all being strongly dominated by emission 
lines. The positions of the most prominent emission fea- 
tures are labelled; i.e., the H-like lines of Cvi, Nvn, Ovm 
and Nex, the strong He-like triplets of Ovil and Neix as 
well as a number of Fe xvn - Fe xix lines. The peak forma- 
tion temperatures of these lines cover temperatures from 
~1.5MK ( Cvi) over -3MK (Ovm) and ~6MK (Nex) 
up to ~ 8 MK (Fe Xix) but plasma with temperatures up 
to 30 MK can contribute to these lines. 

The EPIC PN data, with moderate spectral resolu- 
tion, are displayed in Fig.|3for our target stars; note that 
the EQ Peg EPIC spectra were extracted from an annulus 
because of pile up. In addition to the features mentioned 
in the RGS section the continuum at high energies and 
lines due to highly ionized iron (up to Fexxv at 6.7 keV) 
and H- like/He- like magnesium, silicon and sulfur are de- 
tectable. The observed features cover roughly the temper- 
ature range from 2 to 70 MK in terms of peak formation 
temperatures. Inspection of the PN spectra of the sample 
M dwarfs shows differences in the high energy slopes, the 
obviously shallower slope of the AT Mic spectrum indi- 
cates the higher coronal temperatures of this source. 
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Table 4. Results on temperatures (keV), emission measures (10 cm J ) and abundances relative to solar photosphcric 
values, errors are 90% conf. range, (d.o.f. given in brackets). Lx (10 28 ergs _1 ) indicates the 0.2-10.0keV luminosity. 
The 6-T model grid is 0.2, 0.3, 0.6, 1.2, 2.4, 4.8 keV which corresponds to 2.3,3.5,7.0, 14.0, 28.0, 56.0 MK. 



PN (low) PN (high) | RGS~ 



PN (low) PN (high) 



RGS 



MOS 



PN 



MOS 



PN 



EQ Peg - 3-T model 



EQ Peg - 6-T model 



c 





71 -0 


17 

14 





























70_ 


18 

13 


























IN 


() 


b9_ 


19 
16 





























KS+° 

b8_ 


22 
15 


























{J 





54 -0 


09 
08 





,7+0 
i7 -0 


05 
01 





58 -0 


09 
08 





7fl+° 


18 
14 


() 


AS + ° 

48_ 


10 
07 





59 -0 


11 

06 


-0.05 





58_ 


09 
08 


() 


70+0 


18 
15 





48_ 


08 
07 


Nc 


n 
u 


Qfi + ° 

yb -o 


16 
13 


n 
u 


70+O 

'3_ 


14 
12 


u 


QO+0 

yd -o 


19 
18 


n 
u 


y2_ 


32 
37 


1 


ui_ 


27 
22 


n 
u 


yb -o 


15 
13 


69+ 09 

u-uc7 -0.109 




nn+° 


16 
13 


n 
u 


Qfi+o 
yb -0 


26 
21 




11+° 

1X -0 


20 
17 


Mg 





30_ 


14 
14 





32_ 


08 
07 


f) 


,r + 

25_ 


11 

09 





99+O 
22 -0 


19 
15 





97+° 

27_ 


15 
11 





97 + ° 


15 
13 





33±° 


08 
08 





9S + ° 


11 

09 





22 -o 


17 
14 





00+0 
33_ 


14 
12 


Si 


o 


58+° 


31 
29 





58 +0 


10 

09 


() 


40+° 


12 
10 


o 


55 +0 


23 
18 


o 


29+° 

zy -0 


15 
11 


o 


56+° 


34 
28 





60±° 


10 
11 


o 


45+° 


13 
11 


() 


58+0 


23 
19 





37+° 
°' -0 


15 
13 


S 





27±°o 


18 
15 





68±° 


17 
16 





27+° 


21 
20 





63±g 


10 
34 





ooi° 


18 
00 





28+_° 


2 3 
17 





69±°o 


18 
16 





27+_° 


21 
22 





69+2 


44 
37 





ootg 


17 
00 


Fe 





^t° 


08 
06 





3lt°o 


on 

04 





47t°o 


1)8 
07 





54+° 


16 
12 





43+° 


11 
07 





43t° 


08 
04 





32+0 


05 
05 





A6±° 


118 
06 





5210 


13 
11 





«±s 


09 
08 


kTl 





26±° 


01 
01 





29±° 


03 
02 





26±° 


01 
01 





27t° 


01 

02 





24t° 


02 
02 
































EMI 


1 


G2±° 


32 
28 


2 


81±°o 


81 

53 


1 


8lt°o 


41 
39 


1 


63±g 


81 
48 


2 


06i° 


56 
54 





52 + " 


15 
15 





17±°o 


37 
12 





™±°o 


21 
20 





30±° 


22 
23 





96t° 


37 
36 


kT2 





G3±° 


02 
02 





60±° 


03 
02 





64t°o 
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3.2.2. Spectral fits and coronal temperatures 

Closer inspection of the M-dwarf spectra and their indi- 
vidual lines reveals further differences. To quantify these 
differences the spectra from each type of detector were 
fitted separately with the 3 -T model and the 6 -T model. 
Integration times of the spectra may differ due to vary- 
ing exposure times and the influence of high background 
periods on the individual detectors. We examined the in- 
fluence of flaring on the temperatures, emission measures 
and abundances of the plasma. For this purpose we ex- 
tracted spectra from flaring and quasi-quiescent intervals 
as defined in Tabic [3] and modelled the PN spectra for the 
different phases of activity, keeping all model parameters 
as free parameters. The results of this fitting procedure 
are summarized in Table all errors denote 90 % confi- 
dence range and were calculated by allowing variations of 
all other free model p arameters. Note the c orrected flux for 
EQ Peg compared to lRobrade et all l(2004l) . The quality of 
these fits in terms of the \ 2 test statistic did not improve 



when further components were added to the model. We 
point out that the calculated errors are statistical errors 
only and additional uncertainties are introduced by the 
atomic data and instrumental calibration. 

Comparing the fit results for the different instruments 
one finds good agreement in general, but deviations of 
particular results are obvious. One source of uncertainty 
is based on the interdependence of emission measure (EM) 
and elemental abundances; a higher EM of a temperature 
component leads to the reduction of the abundance of el- 
ements producing the strongest lines in this temperature 
range and vice versa. The difference in fit quality between 
two models is often only marginal. This is especially true 
for EPIC data, where the best-fitting models result some- 
times in higher EM with lower abundances compared to 
RGS data. However, a combination of these properties is 
in good agreement, e.g. a simple scaling of the summed 
EM with the summed abundances of iron and oxygen (see 
Table EJ. Iron and oxygen are the major contributing el- 
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ements of the emission of M dwarfs in the considered en- 
ergy range; their abundances can therefore be determined 
best. Similar results were obtained when using the oxygen 
abundance for the cooler components and iron for the hot- 
ter components. Also the emission measure is in reality a 
continuous distribution which results in an additional in- 
terdependence between temperatures and strength of the 
respective EM components, noticeably in larger, often de- 
pendent, errors affecting neighbouring components. The 
3 -T and the 6 -T model lead to fully-consistent results on 
abundances for all targets; differences are within the errors 
and negligible when compared to the uncertainties arising 
from the use of different detectors. 

One effect of calibration uncertainties is that the best 
fit models of the different detectors do vary in X-ray flux, 
especially below 1.0 keV. Typical RGS models result in 
a ~ 20% lower X-ray flux, mainly in the energy range 
0.2-1.0keV compared to EPIC, whereas the PN models 
are a few percent more luminous than the MOS models. 
The RGS luminosities in Table0]have been extrapolated to 
the 0.2- 10.0 keV band for comparison with the EPIC lu- 



Table 5. Deviations of the XMM instruments from the 
mean scaled summed EM (S £ , M*(24 Fe + J 4 )[10 51 cm -3 ]). 



Instr. 


EQ Peg 


AT Mic 


AD Leo 


EV Lac 


RGS 


-5% 


+3% 


-4% 


-5% 


PN 


+ 10% 


+7% 


+9% 


+7% 


MOS 


-5% 


-11% 


-4% 


-3% 


sc. SEM 


4.51 


15.32 


3.08 


2.97 



minosities. Smaller cross-calibration errors, particularly at 
energies below 0.5 keV, are also seen in observations made 
with the same detector but different setups, especially 
with different filters. Enclosed energy correction from an- 
nular extraction regions, necessary for piled-up sources, 
also influences the measured flux, which can be reduced 
by up to ^10%. Additional effects are expected for the 
X-ray flux of the binaries, but in our case they are dom- 
inated by one component and/or the separation is small 
compared to the used ext raction region . Calib ration efforts 
are still ongoing, see e.g. iKirsch et all $2004) for updated 
XMM-Newton (cross)-calibration. 
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We note that that the EPIC detectors are not very ap- 
propriate to measure changes in the abundance pattern for 
the different phases of activity as they sometimes result in 
unphysical solutions like a zero abundance. On the other 
hand, the RGS is not very sensitive to high temperature 
plasma and in the 6 -T models applied to RGS data some- 
times plasma is put into the hottest components, that is 
not seen in the EPIC models. Since the contributions from 
these components in the RGS are also compatible with a 
zero value, we expect the EPIC models to give a more 
realistic EMD at these energies. Overall we consider the 
abundance determination below ~ 1.5 keV to be more ac- 
curate for RGS data, while properties of the hotter plasma 
are best determined with EPIC data. Therefore we anal- 
ysed additionally for the different phases of activity the 
distribution of elemental abundances with RGS data, and 
the emission measures with PN data using a fixed abun- 
dance pattern. The results of flare-related changes will be 
presented in section l3~3l 

The derived temperature structure is consistent among 
the detectors, independent of the model used, when ac- 
counting for the above mentioned uncertainties. The tem- 
perature distributions of all sample stars peak consis- 
tently around 7 - 8 MK ( 1 ke V = 11.6 MK) with the bulk of 
the plasma being located between 2-20MK for AD Leo, 
EQ Peg and EV Lac, whereas the AT Mic distribution 
extends to somewhat higher temperatures. In the EPIC 
3-T models the high temperature component tends to be 
hotter and/or stronger compared to RGS. 

The EMDs for the total observation of our sample 
stars as derived from the 6-T model with the differ- 
ent detectors onboard XMM-Newton are shown in Fig.0] 
Undetected components are plotted with their upper lim- 
its. The EMDs were scaled with the summed abundance 
of iron and oxygen and are in good agreement when com- 
paring the individual instruments but individual temper- 
ature components are sometimes more pronounced. The 
shape of the distribution looks similar for our sample stars 
but the contributions of the hot plasma differ and are 
strongest for AT Mic, which is also the X-ray brightest 
source in our sample. AT Mic is associated with the 
Pictoris moving group llBarrado v Navascues et all [l999; 
IZuckerman et all 1200 l|L with a probable age of about 10- 
20 Myr whereas the other stars are members of the kine- 
matic young disk population with ages of the order of 
1 Gyr. At this age AT Mic would be a still contracting 
pre-main sequence system, therefore the stars are larger 
with higher Lboi and because they are in the saturated 
regime, they have consequently higher Lx- 

We also investigated joint fits, where data of differ- 
ent instruments, e.g. RGS and MOS, are fitted simultane- 
ously. While combined fits with datasets of different de- 
tectors usually turn out a weighted, intermediate result, 
this is not the case when taking RGS data in the full en- 
ergy range and combining it with MOS data above 0.8 or 
1.0 keV. This choice is reasonable because the relative en- 
ergy resolution of the EPIC instruments improves towards 
higher energies, whereas the effective area of the RGS 



drops above 1.0 keV. The derived abundances in these fits 
increase by ~ 20-30%, whereas emission measure drops, 
compared to the pure RGS model; the derived temper- 
atures change only marginally. The effect on individual 
elements and respectively emission measure components 
slightly depends on the chosen overlap and is similarly ob- 
served with PN data. Discarding the RGS data at higher 
energies, e.g. above 1.5 keV, has only minor influences on 
the derived models for our stars. 



3.2.3. Coronal abundances 

The abundance patterns of our sample stars are shown in 
Fig. where we plot the abundances with respect to solar 
photospheric abundances against the FIP (First Ionization 
Potential) of the corresponding element. Here we show the 
results for the total observation that are constrained best. 
The data of the individual stars are separated by an offset 
of one on the vertical axis and a dotted line is plotted at 
0.5 solar abundance to guide the eye. For elements with 
stronger lines in the RGS we accept these values (i.e., for 
Fe, C, O, N, Ne), for the other elements we use the MOS 
values, which has a better spectral resolution compared 
to PN. The higher oxygen v alues obtained with the solar 
oxygen abundance given by lAllende Prieto et alJ l(200lh 
are also plotted, indicating the current uncertainties in 
the solar abundance pattern. A discussion of problems in 
determining solar photospheric abundances can be found 



determining 
in lHolwegerl 



The determined metallicities of our sample stars vary 
only moderately between the individual stars; we do not 
consider these differences to be very significant. However, 
all abundances are subsolar when compared to solar pho- 
tospheric values. Unfortunately, photospheric abundances 
for M dwarfs are usually unkn own and the rare mea sure- 
ments have large uncertainties. iNaftilan et alJ Ijl9921) pub- 
lished photospheric abundances for AD Leo, arguing for 
a metallicity of ~ 0.4 solar for the analysed low FIP met- 
als; iron was determined to ~ 0.5 solar, in agreement with 
our coronal measurements but chrom ospheric contamina- 
tion could not be ruled out completely. IZboril et alJ ( 1998h 
published a metallicity of ~ 0.6 solar for a number of K and 
M dwarfs from optical/IR observations. These abundance 
values are challenging because our sample stars are mem- 
bers of the kinematic young disk populatio n (with AT Mic 
being even younger), that are fast rotating I Delfosse et all 
1998) and highly X-ray active, which indicates a younger 
age in contrast to the reduced metallicity one would asso- 
ciate with older stars. 

If the photospheric metallicities in our sample stars are 
indeed as low as suggested by the authors cited above, i.e. 
in the range 0.4-0.6 solar, the abundances of low FIP ele- 
ments would actually be consistent with the stellar photo- 
spheric values or possibly slightly depleted, while the high 
FIP elements would actually be enhanced in the coronal 
material. An overall coronal metal abundance deficiency 
would then not be present, rather one would find a situa- 
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tion reminiscent of the solar corona, i.e., an enhancement 
of specific (but different) elements. 

When comparing star-to-star variations of abundance 
as a function of FIP, one finds a tendency that the varia- 
tions of low FIP elements are smaller than those of high 
FIP elements, e.g. neon varies by more than 50 % com- 
pared to iron where the effect is below 10 %. Whereas the 
absolute values differ moderately, the observed abundance 
pattern looks very similar for all stars. In general, abun- 
dance ratios are more robust than absolute values because 
they are nearly independent of the EM. We already men- 
tioned that while low activity stars show the FIP effect 
more active stars have shown none or an inverse FIP ef - 
fect (IFIP) l|Brinkman et all l200li lAudard et all l2003h . 
As a major trend the inverse FIP effect is clearly present 
in all analysed sources but individual elements show occa- 
sional deviations, e.g. carbon, which is mainly determined 
by the Cvi line at the low energy end of the RGS. The 
IFIP effect is detected with all instruments; however, the 
strength of the IFIP or abundances of individual elements 
do vary. This can be mainly attributed to instrumental 
properties, e.g. the often low abundance of sulfur in RGS 
fits. Here the detectable lines lie at the high energy end of 
the detector where the effective area is quite small. The in- 
verse FIP effect is strongest for the most luminous source 
AT Mic and appears to be correlated with activity. Taking 
the abundance ratio Ne/Mg as an indicator of the strength 
of the IFIP effect (since these elements have the largest 
difference in FIP), we find the largest ratio for AT Mic, in- 
termediate values for AD Leo and EQ Peg and the lowest 
ratio is found for EV Lac which also exhibits the lowest 
absolute abundance values. The same trend is observed 
for other ratios like Ne/Fe or O/Fe and likewise in the 
results obtained with the MOS detector. 

3.3. Flare-related changes 

Significant variations in the abundance pattern were not 
detected with the EPIC data. We also searched for vari- 
ations with RGS data, which suggests an enhancement 
of low FIP elements like iron during strong flares, most 
prominent in the large flare on AT Mic when compared 
to the quasi-quiescent phase. While the elements with in- 
termediate FIP show no clear trend, the high FIP ele- 
ment neon may also be somewhat en hanced, a trend al- 
ready found bv iR.aassen et alJ l)2003h . We caution, how- 
ever, that the measurement errors are large and that a 
constant abundance pattern is also acceptable within the 
90 % confidence range derived from the model. This effect 
is still weaker or undetected for the other sample stars. 

In contrast to abundances, emission measures and tem- 
peratures of the coronal plasma change significantly. As an 
example we show the PN spectrum of the AT Mic observa- 
tion in Fig. [HI separated into the three phases of activity. 
Not only does the plasma emission measure increase, but 
the spectra also harden towards higher activity as more 
and more emission is found in the higher temperature com- 



ponents during flaring. This is also seen in the growth of 
the iron lines at 6.7 keV, which are predominantly formed 
at temperatures above 40 MK. 

A quantitative study of the temporal evolution of 
plasma temperature during flaring was carried out via the 
emission measure distribution (EMD). Because the most 
significant changes are expected to be in the high tempera- 
ture components of the plasma, we used PN data to study 
the changes of the EMD. For this purpose we use the 6-T 
model, and keep the abundances fixed at the values as de- 
termined from the total observation to derive the EMDs 
for the different phases of activity. The first two compo- 
nents represent the cool plasma (2-5 MK) , followed by the 
medium temperature plasma (5-20MK) and the last two 
account for the hot plasma (20-70MK). 

The observed sources vary in terms of X-ray luminos- 
ity and its increase during flaring. The largest increase 
is found for the most luminous source AT Mic, where 
the X-ray luminosity during the flare peak phase reaches 
~ 44 * 10 28 ergs -1 ; the source brightened by nearly a fac- 
tor of two for about an hour. AD Leo, EQ Peg and EV Lac 
are about half an order of magnitude X-ray fainter than 
AT Mic. A graphical representation of the results on the 
EMDs of the lower activity stars separated into two dif- 
ferent phases of activity (three for the most active source 
AT Mic) is shown in Fig.0 

The weakest flaring activity is found on EQ Peg, with 
changes in EM affecting mainly the medium temperature 
component. The summed EM of the cool component is 
constant within errors whereas the medium component in- 
creases significantly; typical flare temperatures are around 
10-15MK. The hot component is only marginally de- 
tected, contributing only 3 % to the EM, and does not ap- 
pear to be significantly affected by flaring. An interme- 
diate case is observed for AD Leo and EV Lac; here the 
cool components again stay nearly constant, the medium 
components increase significantly and the hot component 
also increases during flaring. The flare temperatures are 
in the 15-40MK range. Remarkably the existence of two 
more separated parts of the EMDs derived from PN data 
appears also in quasi-quiescence and is similarly found on 
other temperature grids. However, the respective plasma 
contributes only a few percent to the total emission mea- 
sure and the EMDs derived from MOS data are likewise 
more compact in the sense that the plasma appears to be 
more concentrated at medium temperatures. 

The strongest flaring is observed on AT Mic; here the 
cooler plasma below the main peak at ~ 7 MKs is nearly 
unaffected by the flaring. The main increase is due to the 
hottest plasma components and large amounts of addi- 
tional plasma were detected during flaring at temperatures 
between 15 - 60 MK. In the low state a second weaker peak 
is visible in the EMD between 25 and 30 MK, followed by a 
steep decline towards higher temperatures, indicating no 
significant amount of plasma above 40 MK. With inter- 
mediate conditions during the medium activity phase the 
plasma distribution during the large flare exceeds 60 MK. 
The contribution of the hot component to the EMD during 
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the flare is around 25 %. The influences of possible changes 
of the abundance pattern on the EMD are minor and we 
find that the basic behaviour and shape of the EMD is 
comparable for models with free or fixed abundances. 

Our analysis of the EMDs shows that flaring does not 
affect the cool plasma in a significant way. The shape of 
the EMDs are quite alike at cooler temperatures for all 
stars and peak at temperatures around 7-8MK whereas 
the absolute value correlates with X-ray luminosity. The 
main influence during flares is found on the amount and 
temperature of the hotter plasma: stronger flaring and 
more active stars produce more and hotter plasma. As 
to the emission of the stars in the quasi-quiescent state 
a relation between X-ray luminosity and average plasma 
temperature is suggested while the overall shape of the 
EMDs is again similar. Only the distribution of binary 
system EQ Peg, whose E MD is created by two active and 
slightly different stars (see iRobrade et all l|2004(l 1. appears 
broader in all detectors. A hot tail is quite common, some- 
times with a second smaller peak, extending to ~ 30 MK 
in the EMD. This may indicate the existence of sepa- 
rated plasma components at different temperatures which 
are located in different coronal structures even in quasi- 
quiescence, i.e. the main part of the plasma in cooler, more 
quiescent structures and a smaller amount in hotter and 
more active structures. These structures would also be a 
natural explanation for the constant flickering observed in 
the light curves. 

3.4. Loops 

Assuming quasi-static loop-like structures as building 
blocks for the coronae of M dwarfs the derived EM dis- 
tributions suggest coronae composed out of different loop 
components with different temp eratures. Using the RTV 
scaling law l)Rosner et al.Lfl978l) for simple loops in hydro- 
static equilibrium and with consta nt pressure and c ross- 
section, and the densities derived bv lNess et al.l (|200J), we 
are able to calculate the size of the loops. In cgs units we 
obtain: 

with the electron density n e , the loop top temperature 
T m ax and the (semi) loop length L. The typical electron 
densities of our sample stars in log(rt e ) as measured from 
the Ovn triplet are around 10.5, and for the Neix triplet 
in the range 11.1-11.5, but for Neix the errors are about 
0.5 translating into a factor of three in linear units. The 
Ovil triplet traces plasma around 2MK, the Neix triplet 
plasma around 4MK, thus indicating denser plasma at 
higher temperatures. However both triplets trace only the 
relatively cool plasma since we find average coronal tem- 
peratures of about 10 MK for the less active stars and 
15 MK for AT Mic from the EMDs, in good agreement 
with the values calculated with a scaling law correlat- 
ing X-ray lu minosity and stella r radius (1?*) with average 
temperature l|Ness et al.Ll2004^l . While the assumptions of 
quasi-static loops and the measured densities can be well 



justified for the cooler and the quiescent components, re- 
sults for the hotter plasma should be taken as upper limits, 
especially when taking into account the plausible higher 
density of the hotter plasma. 

Assuming stellar radii of 0.35/0.3 solar radius for spec- 
tral type M3.5/M4.5 and for the binaries for simplicity 
one dominant component, we determine loop lengths for 
the coolest plasma (2-5MK) below 0.01 1?* on all stars. 
Taking the mean plasma density derived from neon, the 
loop lengths are about 0.03-0.071?* for the average coro- 
nal temperature of our stars a nd increase up to .3 1?* for 
the hotter plasma components. iNess et all l|2004|) find the 
density of the hotter plasma consistent with the low den- 
sity limit (log 716^12.7) for the used line ratios of FeXXl, 
so an higher density as the one derived from Ne IX in the 
hotter loops would further scale down the determined loop 
length according to the above scaling law; instead of an in- 
creased loop size, an increased density also can explain the 
observed temperatures. Repeating the calculations with 
the low density limit of Fexxi, we derive lower limits 
for the length of the hotter loops, which are scaled down 
by a factor of 20 compared to the above derived results. 
Therefore we conclude that the dominating coronal struc- 
tures for the cool plasma are very small, i.e. L< 0.011?*, 
for the medium temperature plasma small in the sense 
that L < 0.11?* and still co mpact (L^0.31?*) for the hot- 
ter plasma observed here. iFavata et alJ (|2000al) studied 
AD Leo data from various missions, analysed decay times 
of flare loops and likewise found its corona to be compact, 
even for flare structures. 

In this picture the cool emission is produced in rel- 
atively small loops nearly unaffected by flaring, while 
medium and hot temperature loops are larger, probably 
have higher densities and are more numerous on more ac- 
tive stars or during flares. 

3.5. Comparison with other results 

We compare our results with other measurements and 
earlier works to check the consistency of the derived re- 
sults and to investigate long-term behaviour. We find that 
the results presented here agree well with other analyses, 
which often used different techniques or dealt with indi- 
vidual objects only, leading to a rather consistent picture 
of the X-ray properties of active M dwarfs as a class and 
demonstrate the robustness of the results. 

Comparison of the measured fluxes of o ur sample stars 
with e.g. the ROSAT all-sky survey values ijHunsch et all 
1999J) shows agreement to within a factor less than two. 
This is of the same order as the variations found dur- 
ing the moderate activity in the XMM-Newton observa- 
tions, demonstrating the stability of the X-ray emission. 
Favata et alJ l)2000a|) studied the quasi-quiescent corona of 
AD Leo from various missions and found the X-ray lumi- 
nosity to be remarkable constant with variations less than 
a factor of two, suggesting the absence of strong cyclic 
variations; our results agree well, expanding the time pe- 
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riod of the observations to over 20 years. The X-ray lumi- 
nosity of EQ Peg determined in this work agrees likewise 
with the measurements made with Einstein more than 20 
years ago. Combined, these findings establish long term 
stability of the quasi-quiescent corona for active M dwarfs 
in general. 

Comparing our results with previous analysis of the 
XMM-Newton data, our AD Leo results are in good agree- 
ment with those obtained bv lBesselaar et akl 1 200.*^ . who 
also applied the global fitting approach, especially when 
considering the different spectral code and model, atomic 
database and of course the X-ray dataset used in the 
analyis. The same is tr u e for the AT Mic results pre- 
sented in iRaassen et al.l l|2003j) with a similar analysis; 
note that errors stated in these publications are only la. 
Discrepancies concerning the absolute values of emission 
measure and abundances are explained by their mutual 
interdependence, and a combination, e.g. the product, of 
these properties is in much better agreement. Concerning 
the shape of the DEM, which is determined there with 
polynomial methods from MOS and RGS data, results are 
comparable with our EMD based on multi-temperature 
fits of PN data. 

iMaggio et al.l l|2004f) presented a line-based analysis of 
AD Leo data taken with Chandra LETGS and also found 
the inverse FIP effect. However, they determined the iron 
abundance from the line-to-continuum ratio in the low- 
count spectrum to be almost as high as that in the so- 
lar photosphere. Consequently, nearly all other elements, 
which are scaled with the iron abundance, appear t o be 
supersolar. On the other hand, iFavata et al ] ll20004 de- 
termined the iron abundance from ASCA data to ~0.2 
solar, low but comparable to our EPIC values. While 
the abundances determined with medium resolution spec- 
troscopy are apparently often low, the results obtained 
with the global fits of hig h resolution data in our work 
or in lBesselaar et all (|2003l) appear more reliable. We also 
point out that our joi nt fit approach which uses a data se- 
lection comparable to lAudard et al.l (|2004^ . who obtained 
consistent results on the abundances with line-based and 
global models applied to high resolution X-ray spectra, 
results in slightly increased but still significant subsolar 
abundance values. 

Whereas the global metallicitics were not determined 
uniquely in the various analysis, the inverse FIP effect 
is consistently found in all analysis of active M dwarfs, 
despite the fact that different spectral models and atomic 
codes were applied to medium and high resolution X-ray 
data, demonstrating that the inverse FIP effect is certainly 
a robust result. 

4. Summary and Conclusions 

We have carried out a comparative analysis of XMM- 
Newton observations of M dwarfs with spectral type 
M3.5-M4.5, i.e., AT Mic, EQ Peg, AD Leo and EV Lac. 
Light curves revealed frequent flaring of all observed 
sources on different timescales and with various strengths 



of the individual flares. While AD Leo, EQ Peg and 
EV Lac are comparable with respect to activity and X-ray 
brightness, the probably younger system AT Mic is half 
an order of magnitude X-ray brighter. From X-ray spec- 
troscopy, abundance patterns and emission measure dis- 
tributions were derived for the sample stars and for differ- 
ent states of activity. We find that many X-ray quantities, 
like the coronal temperature structure and the abundance 
pattern of the analysed M dwarfs are very similar despite 
the fact that the X-ray luminosity and flaring behaviour 
differs significantly. 

All analysed M dwarf coronae show the inverse FIP 
effect, demonstrated in this work for the first time for 
EQ Peg and EV Lac, establishing this effect for the coro- 
nae of active M dwarfs as a class. The determined abun- 
dances are lower than solar photospheric abundances; 
however, the few available measurements of stellar pho- 
tospheric metallicities for M dwarfs also point to subso- 
lar values, comparable with our coronal measurements. 
This suggests that a strong coronal metal abundance de- 
ficiency may actually not be present in our sample stars. 
Properties like X-ray brightness, temperature of the flare 
plasma and thus hardness of the X-ray spectrum are high- 
est for the most active source AT Mic, and increase during 
flares for all observed sources, suggesting a correlation be- 
tween these properties and activity. The strength of the 
inverse FIP effect is also greatest in AT Mic, but does not 
increase during flaring. Time-resolved abundance analysis 
suggests that preferably low FIP elements are enhanced 
during flares, but this finding is not statistically signifi- 
cant. 

We find from the derived EMDs of the stars that the 
plasma temperature structures of the coronae are sim- 
ilar, with a maximum at 7-8 MK and a hot tail with 
a stronger contribution of hotter plasma in the X-ray 
brighter stars. The emission measures of the hotter plasma 
components increase significantly during flares and higher 
temperatures are reached in stronger flares, while the cool 
plasma is only marginally affected by flaring. Assuming 
that the coronae are composed of loop-like structures, this 
indicates a distribution of structures in which the cooler 
plasma is confined in many smaller loops while other, pre- 
sumably larger and denser loops, dominate the additional 
hot emission generated by flare events. The coronae are 
found to be compact, i.e. the largest loop lengths are 
smaller than one third stellar radius. 

From comparison with older measurements we con- 
clude that the global X-ray properties show long-term sta- 
bility while short time-scale behaviour is characterised by 
frequent variations on top of a level of 'baseline activ- 
ity'. With the XMM-Newton data we are able to investi- 
gate the X-ray properties of stellar coronae to a previously 
unknown precision. As these observations provide only a 
snapshot in the life of a few M dwarfs, larger samples or 
longer observations may show more extreme or different 
phenomena. Nevertheless, the consistency and congruence 
makes our results appear typical for active M dwarfs. 
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Fig. 2. RGS spectra of the sample stars, the most promi- 
nent emission features are labelled. 
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Fig. 3. PN spectra of the analysed M dwarfs; top to 
bottom: AT Mic (red), AD Leo (green), EV Lac (blue), 
EQ Peg (black). [Colour figure in electronic version.] 



Fig. 1. Light curves of the binaries EQ Peg and 
AT Mic and the single stars AD Leo and EV Lac; PN 
(+MOS2/red) data, 100 s binning. Frequent flaring is ob- 
served on various timescales; selected flaring periods are 
indicated by horizontal arrows. 
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Fig. 4. EMDs as derived with the 6-T model and scaled 
with the summed abundance of iron and oxygen, total 
observation each (RGS: Asterisks, MOS: Diamonds, PN: 
Squares). 



A AT Mic 
□ EQ Peg 
O AD Leo 
'X. EV Lac 



L_MgFeSi 



$ § 
S C 



■f-r 

N 



s 

Ke 



FIP (eV) 

Fig. 5. Elemental abundances in stellar coronae relative 
to solar photospheric values vs. FIP for our sample stars; 
an inverse FIP effect is present in all stars. Each star is 
shifted by one solar abundance on the vertical axis, the 
dashed lines indicate 0.5 solar abundance in each case. 
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Fig. 6. Spectra of AT Mic during different states of ac- 
tivity, PN data from selected time-intervals as defined in 
Table |3 (high: black, medium: green, qq: red). [Colour fig- 
ure in electronic version.] 
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Fig. 7. EMDs of our sample M dwarfs in different states 
of activity derived from PN data. The upper plot shows 
the flaring phase, for AT Mic the medium activity phase is 
additionally plotted as a dashed line while the lower plot 
shows the quasi-quiescent state. Undetected components 
are plotted with their upper limits, denoted by downward 
arrows. 



